Frequency Evolution of Neutron Peaks Below T^. Commensurate and Incommensurate Structure in 

LaSrCuO and YBaCuO 



O 
O 

o 

(N 
X) 



o 
o 

u 

Oh' 



X3 
O 

o 



> 

O 

m 

00 

o 

On 
On 
-I— > 



I 

O 

o 



X 



Ying-Jer Kao 

James Franck Institute, University of Chicago, Chicago, Illinois 60637, USA 

Qimiao Si 

Department of Physics, Rice University, Houston, TX 77251, USA 
K. Levin 

James Franck Institute, University of Chicago, Chicago, Illinois 60637, USA 
(February 7, 2008) 

We study the evolution of the neutron cross-section with variable frequency uj and fixed T below Tc in two dif- 
ferent cuprate families. Our calculations, which predominantly probe the role of d-wave pairing, lead to generic 
features, independent of Fermi surface shapes. Among our findings, reasonably consistent with experiment, are 
(i) for Lo near the gap energy A, both optimal LaSrCuO and slightly underdoped YBCO exhibit (comparably) 
incommensurate peaks (ii) peak sharpening below Tc is seen in LaSrCuO, (iii) quite generically, a frequency 
evolution from incommensurate to commensurate and then back to incommensurate structure is found with in- 
creasing UJ. Due to their naiTow u regime of stability, commensurate peaks in LaSrCuO should be extremely 
difficult to observe. 

PAC numbers: 74.25.Ha, 74.25.-q, 74.20.-z 



The field of neutron scattering in high temperature super- 
conductors has seen a variety of recent experimental discover- 
ies associated with incommensurate and commensurate struc- 
ture at r < Tc These neutron data remain of central 
importance in the field of high Tc superconductivity: the com- 
bined momentum and frequency range covered by this tech- 
nique is wider than that in virtually all other spectroscopies. 
It is mainly through this extended range that these neutron 
results have direct bearings on such important issues as dy- 
namical stripes and where the condensation energy comes 
from §1. 

The goal of the present paper is to systematically address 
these observations over the entire momentum and frequency 
range that the experiments have covered, in the two different 
cuprate families (LaSrCuO and YBaCuO). In the process, we 
show that all of the above commensurate and incommensu- 
rate features are compatible with d-wave pairing superposed 
onto the normal state Fermi surfaces. In contrast to T > Tc, 
the details of the Fermi surface shape are of relatively less 
importance, and serve primarily to select out the oj regime 
where various commensurate or incommensurate features can 
be observed. The importance of the present work derives from 
the panoply of different experimental observations which are 
semi-quantitatively addressed here. These calculations have 
no adjustable parameters(besides those which were used orig- 
inally 10 1 to fit some aspects of normal state data), so that 
their success or failure, upon comparison with experiment, 
should help select out viable theories of the cuprates. 

Here, following previous work above ^ and below [ |ic| , [ri| ] 
Tc, we apply an RPA, three band scheme which we have de- 
veloped to treat the effects of very strong Coulomb correla- 
tions [p^. It should be stressed that our RPA approach is 
not a weak coupling RPA. The Lindhard function Xo used 
here is appropriate whenever the spin excitations are asso- 



ciated with spin 1/2 and with an underlying Fermi surface- 
be it in a Fermi liquid |pj[-p^], or in a spin-charge separated 
metal [p^|-po||. Our starting point is the dynamical suscepti- 
bility [0 x(g,w) = x°iq,^)/['^ + J{q)x°{q.^)] where, at 
low T, the dominant contribution to the imaginary part of x° 
is given by 



Here u{k,q) = (1 - {£,kik+q + ^k^k+q) / EkEk+q), and 
represents the "bare" particle energy relative to the chem- 
ical potential while Ek is that of the superconducting quasi- 
particles. The important function E2{k, q) = Ek + E^+q will 
play a key role in our analysis. 

The RPA neutron cross section reflects a competition be- 
tween effects associated with the Fermi surface shapes and 
pairing symmetry [via Xo] and those from the residual ex- 
change interaction [ p^ J{q) ~ Jo[cos qx + cos qy], which de- 
rives from Cu-Cu interactions via the mediating oxygen band. 
While the YBaCuO system is a two layer material, our past 
experience has shown that most of the peak structures 

associated with the neutron cross section are captured by an 
effective one layer band calculation, which we will investi- 
gate here. For definiteness, we have fixed the temperature at 4 
K and assume the electronic excitation gap to be described by 
an ideal rf-wave, A{q) = A(cos q^ + cos qy), where at T = 4 
K, A is taken to be 17 meV in YBaCuOe.e and 8 meV in op- 
timally doped LaSrCuO. The breakdown of the Fermi liquid 
state is addressed only insofar as there may be precursor pair- 
ing orpseudogap effects, which lead to an excitation gap in the 
Lindhard function Xo above Tc. Our calculations were based 
on a numerical procedure in which the Brillouin zone is sub- 
divided into tetrahedral microzones [pH|. Our peak heights 
are represented in arbitrary units, which are best quantified by 
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FIG. 1. Frequency evolution of neutron cross section (Imx) for 
YBaCuOs.e- Units are arbitrary witli corresponding normal state val- 
ues given as a basis in text. Lower panel indicates the position of the 
dominant peaks in reciprical-lattice units (r.l.u.). 



noting that in the normal state, the peak values are around 2 for 
LaSrCuO and for underdoped YBaCuO, whereas for the latter 
compound at optimal doping, the value is around 1 . Presum- 
ably all peak structure with intensity less than this is not cur- 
rently observable (since the normal state peak of YBaCuOy is 
essentially absent experimentally [p2[)). 

In Figure 1 we show the evolution of the neutron peaks 
in underdoped YBaCuO as a function of uj. These incom- 
mensurate peaks are first seen at « A; as frequency in- 
creases, the incommensurability is found to continuously de- 
crease. This decrease is most apparent in the immediate vicin- 
ity of the onset of the resonance, [or (tt, it) peak] which can 
be read off from the lower part of Figure 1 , to be at around 
27 meV, somewhat less than 2A. Just above resonance, the 
(tt, tt) peak becomes flat-topped possibly weakly incommen- 
surate. It then broadens and remains structureless between 30- 
40 meV. Finally, above 45 meV, clear incommensurate struc- 
ture re-appears. We find that our low ui incommensurate peak 
heights are in the ratio of about 1:10 when compared with 
the resonance feature. However, when the integrated spectral 
weight is considered, the ratio is about 1:2. Experimentally 
[jl|], the ratio of spectral weights is found to be 1: 3.8 . 

The lower panel shows the detailed frequency evolution of 
the dominant peak position quantified as (vr, (1 ± 6)tt) ~ 
(0.5, 0.5 ± 5/2) = (0.5, Q) in reciprocal-lattice units (r.l.u.). 
The peaks evolve much as is seen experimentally [^^. The 
primary difference between our observations and these partic- 
ular experiments [^] is that over a range of frequencies, the 
incommensurate features coexist (although, not explicitly in- 
dicated in the lower panel) with the more dominant resonant 
peak. By contrast, experimentally, an energy scale Ec is asso- 
ciated with the frequency at which the various peaks merge. 
It should be stressed, however, that here we have not incorpo- 
rated resolution limiting effects which may affect this detailed 
comparison between theory and experiment. A very early pre- 
diction for this lowest energy scale incommensurability was 
presented by our group in Ref. [p^, where it was shown to be 
a consequence of d-wave pairing and relatively independent 
of the fermiology. Subsequent insights, using a related but 
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FIG. 2. Frequency dependence of ImXo and ReXo for 
Q — (tt, tt). Upper inset is comparable plot for Q = ((1 + 5)n, tt). 
The key characteristic energy scales (see text) are indicated. Lower 
(left) inset shows E2(k, q — Qo) ; here A corresponds to the mini- 
mum of E2{k, q), (ilo) and B to the saddle point of E2{k, q), (oj^h)- 

differently motivated formalism, were provided in Ref [p^], 
which explicitly showed the influence of d-wave pairing on 
ImXo- 

To understand the origin of the various commensurate and 
incommensurate structures, in the main body of Figure 2 we 
plot the prototypical behavior of ReXo and IniXo for q at the 
commensurate point Qo = (tt. tt); in the inset is shown the 
analogous plot for the incommensurate peaks. These plots, 
which were chosen to coiTespond to YBaCuO, contain, in ef- 
fect, a summary of the key energy scales | p3| ] which appear in 
Imxo- Here we emphasize how they are reflected in ReXo- 
These ReXo effects are essential because, through the RPA 
denominator they serve to greatly enhance a given character- 
istic feature in ImXo- The four important energy scales which 
determine the behavior of Imx (via simultaneous effects on 
Imxo and ReXo ) are given as (i) ujo{<i) = min E2{k, q), (ii) 
ujyh, the saddle point of E2{k, Qo) (see point B in the inset), 
(iii) oJcoh the frequency where the rate of change of u{k, q) 
drops abruptly, and, (iv) the onset for commensurate peaks 
ilo = ij-'o(Qo) (see point A in the inset). 

The related implications for ReXo are illustrated in Figure 
2. The onset frequency flo is accompanied by a substantial 
growth in ReXo- However, once the frequency reaches the 
two-particle Van Hove energy ujyh, determined by the saddle 
point shown as B in the inset, ReXo shuts down, as does the 
resonance. Figure 2, thus, shows that because of these ReXo 
amplifications along with the g-structure of J{q), the com- 
mensurate peak will tend to dominate incommensurate struc- 
ture in the range flo < uj < ojvh, because of the small size 
of the RPA susceptibility which, in turn, yields a pole-like 
behavior in Imx- Incommensurate peaks only appear at the 
lower frequencies because their threshold ujo is less than that 
of the commensurate structure. The appearance of incommen- 
surability is associated with the d-wave state and is best seen 
pictorially in the upper right inset of Figure 3, shown for the 
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FIG. 3. Origin of peak sharpening in LaSrCuO below Tc. Super- 
conducting and normal state processes contributing to cross section 
at indicated point are shown in right and left hand figures of upper 
panel, with line intensity corresponding to associated weight. Lower 
panels show the corresponding cross sections. 



LaSrCuO family which contrasts nesting processes in the nor- 
mal state (left) and d-wave superconducting state (right). Here 
the spectral weight of a given process is indicated by the in- 
tensity of the various lines. The magnitude of the incommen- 
surability of the peak is determined by (i) energy conservation 
through the delta function in Eq. (1) along with (ii) coherence 
factor [u{k, q)] effects, which select out the most favorable 
regions from otherwise equivalent nesting vectors. 

To summarize, we can consolidate our observations, along 
with a collection of some of the mechanisms 1 16 , rT| , p4| ] which 
have been proposed for the resonance in YBaCuO into a sin- 
gle inequality: low uj incommensurate peaks appear for ujo < 
ijj < flo, and commensurate peaks appear for flo < uj < ujvh- 
At still higher frequencies the peaks are again incommensu- 
rate, flo is always somewhat less than 2 A because of the nodal 
structure of the cZ-wave gap. 

Indeed, the results of Figure 2 are rather generic and can be 
used to address the LaSrCuO family as well, where unusual 
incommensurate structures in the cross section have been ex- 
perimentally reported below Tc for a range of low uj, above 
A. These, and related "spin gap" features have not yet been 
addressed theoretically. Here, we demonstrate that, in con- 
trast to Ref. [^, these features- which are closely connected 
to the (low energy scale) incommensurate peaks discussed 
above for YBaCuO- are unrelated to fermiology effects and 
to presumed "incommensurate spin structure". They depend 
exclusively on the d-wave pairing symmetry. Figure 3 shows 
a comparison of the cross section in the normal ( left) and su- 
perconducting (right) phases of optimally doped LaSrCuO at 
low w. The peak sharpening below Tc was first observed ex- 
perimentally [0]. Its origin can be traced to the differences 
(above and below Tc) in the initial and final scattering pro- 
cesses along constant energy contours, which contribute to 
Imx" at a given point. These processes are shown in the up- 
per panels of Figure 3 for a particular point q indicated on 
the normal state cross section. It can be seen that in the su- 
perconducting phase, as a consequence of the opening of the 
gap, there is relatively little weight at the q point in question. 
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FIG. 4. Imx for LaSrCuO in the superconducting state along 
q = (0.56^, 0.44/1) (r./.M.) at ft = 1,1.14,1.2. Lefthand inset 
shows the "spin gap" obtained by fitting Imx to a Lorentzian. Right- 
hand inset compares Rex at = in the normal and superconduct- 
ing states. 



thereby leading to a greatly reduced scattering in the region 
between the incommensurate peaks below Tc- 

In order to gain more experimental insight into the behav- 
ior of the neutron peaks of LaSrCuO in the superconduct- 
ing phase, the authors of Ref. [Q| measured Rex at low fre- 
quencies and at the incommensurate points. They also ad- 
dressed the details in the frequency onset of the cross sec- 
tion, measured via Imx, for a range of different wave-vectors. 
To determine Rex, the cross section was fitted to a simple 
Lorentzian-like form [Q] peaked at Ag Here, we use the di- 
rectly calculated Rex to extract Ag. The resulting experi- 
mental curves look qualitatively similar to the theoretical plots 
shown in the right hand inset of Figure 4. A reduction in 
Rex{oj — 0) between the normal and superconducting phase 
can be associated (via general Kramers Kronig relations) with 
a low UJ suppression in the spectral weight of Imx- In the 
present theoretical scheme, this suppression derives from the 
opening of a c?-wave superconducting gap. It is important, 
however, to stress that the d-wave symmetry in Imx is not 
immediately evident from either the plots in Figure 4 or the 
data. Indeed, the main portion of the figure shows Imx ^'^ ^ 
for three different values of the incommensurate wave-vector 
q = (1.12/i7r,0.88/i7r) = (0.56/i, 0.44ft.) (r.l.u.), and indi- 
cates that the onset frequency in Imx is relatively constant in 
q. This onset frequency or " spin gap" A^, is plotted more 
completely for a range of wave-vectors in the upper left in- 
set. Here the range is taken to coincide with its experimental 
counterpart and the results shown compare favorably with ex- 
perimental data from Ref. While the calculated spin gap 
magnitude is larger by a factor of 2 than experiment, it should 
be stressed that the downturn at large h is seen experimentally, 
and is predicted to be systematic. 

Why then is there no d-wave signature in the spin gap fre- 
quency? This follows because the q structure associated with 
the d-wave gap appears inside the integral over k in the Lind- 
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hard function Xo- Moreover, the q points chosen in the main 
body of the figure and in the left hand inset, while consistent 
with those measured in Ref. do not reflect most directly 
the superconducting gap. In our calculations the spin gap As 
is found to be associated with the coherence (coh) factor en- 
ergy scale, oJcoh, which is relatively wave-vector independent 
in the region studied. However, it is the slightly lower, thresh- 
old frequency at which Imx first becomes non-zero ( U0o{q)), 
which most closely reflects A^. (Indeed, lOo varies by a factor 
of 2 for the q values indicated). By contrast with ojcoh, at uoo, 
the cross section has relatively little detectable weight. 

What, then, are the similarities between the LaSrCuO and 
YBaCuO families? There are claims in the literature that 
the incommensurabilities for near-optimal LaSrCuO and for 
YBaCuOg.e are the same, as are the spectral weights once u 
is scaled by the appropriate T^. To address these more quanti- 
tative issues we first note that the w = 23 meV incommensu- 
rability shown in underdoped in YBaCuO in Figure 1 is com- 
parable, w 0.1 (rl.u.), to that plotted in Figure 3 for optimal 
LaSrCuO. 

Moreover, as an estimate of the q-integrated spectral 
weight, we find that the Brillouin zone averages, called 
{Imx{(l,uj)) BZ are equal for YBaCuOg.e and optimal 
LaSrCuO at 22 meV and 13 meV, respectively. If we rescale uj 
by Tc in each case (T^ ~ 60K for YBaCuOe.e [Hi, 35 K for 
Lao.85Sro.i5Cu04 [p7||), it follows that at about the same re- 
scaled frequencies, these intensities are nearly equal. It should 
be stressed that we regard the underlying physical origin of the 
incommensurate peaks in the superconducting state of the two 
cuprates to be the same (as discussed in the context of Figure 
2), but that the characteristic energy and wave-vector scales 
need not be precisely equal. Indeed, the incommensurability 
is now found to be uj dependent [^], which makes the numeri- 
cal comparisons somewhat less meaningful. Finally, it should 
be noted that the main body of Figure 4 is rather generic (as 
seen in Figure 2). These plots of Imx together with their nor- 
mal state counterparts, are in reasonable agreement with the 
results of Ref. 

What then is the most significant difference between the 
LaSrCuO and YBaCuO families? We find that this differ- 
ence derives from the fine details of the bandstructure which 
are then reflected in the range of stability of commensurate 
peaks. The former compound is presumed to have a signifi- 
cantly smaller next nearest neighbor hop ping in tegral it') so 
that the two particle Van Hove feature [|l6,ll| uji,h appears 
much closer to Vlo. Indeed, for LaSrCuO essentially all energy 
scales (for commensurate structure) are compressed down to 
Q,o- In this way there is a much smaller frequency window 
(« 2meV) for commensurate peaks in the LaSrCuO family 
than in its YBaCuO counterpart. Thus far, searches within 
the broader 5 meV interval have failed to see them. While 
this narrow regime of stability will make them difficult to see, 
commensurate peaks appear to be a fairly general feature of 
the RPA- like d-wave approach to the neutron cross section. 

In summary, it should be stressed that the frequency evo- 
lution from incommensurate to commensurate and then to in- 
commensurate peaks in the neutron cross section (as shown by 
the lower panel in Figure 1) is a generic feature of the present 



approach, and because it appears to be observed experimen- 
tally in at least two cuprate families [j!)^, it will be important 
to establish whether an alternative scheme, such as the "stripe" 
picture can lead to similar behavior. 
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